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Optical measuring system 



The present invention relates to a method of measuring amplitude and phase 
variations within a beam of li^t, to a metiiod of characterising a beam of light 
to apparatus for measuring amplitude and phase variations in a beam of light, 
and to apparatus for characterising a beam of light. 

In many optical systems throughput efficiency, crosstalk and noise level 
depend on the alignment of the various components with respect to one another 
and to the optical beams passing through. An important problem during 
assembly of niany free space optical systems is that the beam is not visible and 
its properties cannot be measured as it is incident on the various components in 
the system, except where the component happens to be a detector array. Even ■ 
with a detector array measurements are restricted to the intensity when it is the 
phase distribution that would be more usefully measured. 

Usually measurements are confined to overall coupling efficiency into the 
intended output or some intermediate output. 

The more 'components there are to adjust the more variables there are to 
optimise, and working **blind'' the whole procedure can be very time 
consuming and expensive. A furtiier problem is that the beams in a real optical 
system, although maybe Gaussian-like in theory, often have sidelobes due to 
aberrations, especially when lenses are misaligned. A search for a global 
optimum can lead to an apparent optimum alignment that instead results in a 
sidelobe being coupled into the output, rather than the main peak. The 
existence of such subsidiaiy maxima makes it more difficult to optimise the 
alignment. 




Often liie components themselves are not ideal. For example lenses have a 
tolerance in dimensions and focal length that makes pick and place assembly 
inappropriate, while angle-poMshed fibres have a tolerance in the polish angle. 
As a result the relative longitudinal spacing, transverse offset and tilt of an 
optical fibre and lens may need to be adjusted to suit the properties of that lens 
and fibre. The relative orientation of tiie lens-fibre assembly may need to be 
adapted to the rest of the system. 

Liquid Crystal oyer SiUcon spatial Ught modulators (LCOS SLM) are 
pixeUated devices which may be used for applying phase modulation to 
incident light beams. LCOS SLMs may be used to cany out many optical 
processing functions such as correlation, monitoring a beam by tapping off a 
smaU fraction of the incident energy, routing a beam, changing a beam focus, 
aberration correction, changing a beam shape or changing tiie power carried by 
a beam • 



Bnbodiments of the invention are aimed at providing an abiUty to assist with 
the ahgnment. problems mentioned above by enabling a more complete 
characterisation of beams m an optical system. The invention is however not 
restricted to this area. 

According to one aspect of the invention there is provided a metiiod of 
measuring ampUmde and phase variations in a beam of hght comprising 
causing tiie beam to be incident upon a spatial array displaying a pixellated first 
phase distiibution, in a measuring region of said spatial airay causing the phase 
distiibution to assume a new value wMle retaining the first phase distribution 
outside the measuring region, in the Fourier plane determimng the change in 
intensity resulting from tiie change in phase distribution. 
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According to another aspect of the invention fliere is provided a method of 
characterising a beam of light, comprising disposing a LCOS SLM in the path 
of the beam; causing the LGOS SLM to display a first hologram pattern; at a 
location in said beam where the amplitude and phase of the beam are to be 
5 characterised, changing the hologram pattern to a second hologram pattem; and 
measuring the effect of said change by measuring an intensity 

In an embodiment, the output from the SLM is measured in the Fourier plane to 
detect the Fourier output. 

10 

In an embodiment, the method comprises measuring the intensity in a region of 
the Fourier plane where Fo(x,y) is very weak but g(x,y) is relatively stronger, 
varying the position on the SLM where the perturbation is apphed to form a set 
of measurements of (g(x,y}f(uo,Vo)f taking the square root of , these 
15 measurements to derive values for the relative field amplitude at these 
positions. 

In an embodiment, the method comprises stepping through a sequence of phase 
distributions. 

20 

In an embodiment, the method comprises varying the phase shift in a respective 
siagle pixel. 

According to another aspect of the invention there is provided apparatus for 
25 naeasuring amplitude and phase variations in a beam of Ught, the apparatus 
comprising a pixellated spatial array, each pixel being controllable to apply any 
of plural phase shifts to input light, whereby the array displays a desired 
distribution of phase modulation, means for causing the array to display a first 
selected distribution of phase modulation; means for changing the first 
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distribution in a measuring region of said spatial array to assume a new 
distribution while retaining the jHrst phase distribution outside the measuring 
region^ means disposed in the Fovtrier plane for determiDing a change in 
intensity of li^t resulting from the change in phase distribution. 

According to a further aspect of the invention there is provided apparatus for 
characterising a beam of Kght, comprising a LCOS SLM arranged so that a said 
beam of light can be incident upon it; means for causing the. LCOS SLM to 
display a first hologram pattern; means for changing the hologram pattem to a 
second hologram pattem at a location in said beam where the amplitude and 
phase of the beam are to be characterised; and means for measuring an 
intensity of ligjit to detemaine the effect of said change of hologram pattem. 

In an embodiment, the means for measuring is disposed in the Fourier plane to 
detect the Fourier output. 

In an embodiment, the apparatus further comprises a lens for .providing the 
Fourier ou^ut 

In an embodiment, the apparatus further coinprises a nodrror for providing the 
Fourier output 

There follows a description of how one or more SLMs may be used to aid the 
assembly of the optical system in which they are to act a^ an optical processing 
device, with reference to Figure 1, which shows an exemplary set up of a 
system for use with the invention. 



Referring to Figure 1, an arrangement is shown with a 1-D reflective phase- 
modulating PLC SLM with 540 pixels. The arrangement is a 2-f system with a 
Fourier lens^ and furlh©: comprises a siJica-on-silicon waveguide array. 

A 2-D SLM may be used instead of a 1-D SLM. Use of a reflective SLM is not 
fundamental to the invention. 



Let (u,v) be the co-ordinate system at the SLM and let the incident beam be 
f[u,v) exp i ^M,v> where ^w,v^ describes the amplitude and^u,v) describes the 
phase. As an example the incident field could be an off-axis normally incident 
defocused Gaussian beam. For this example, the incident beam may be 
described by equation (1). 



/(«, V) = exp^ l^'^-^^cMv-v^J-j 



cxpi^M = exp~zfc| '"^^^^ ^^2lt ''''''' I 



where uij^c and Vjj^c the co-ordinates at the centre of the Gaussian beam and 
R is the radius of curvature. 

Consider the SLM to apply a known hologram pattern Ho(u,v) to the iacident 
beam where in general Hq will be a complex function describing phase and/or 
amplitude modulation. In the general case the output field from the SLM is 

Ho(u,v)f(u,v) exp i 0(u,v), 



Measure the output from the SLM in the Fourier plane, using a suitably 
positioned lens or lenses and one or more optical receiving devices to detect the 
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Fourier output at one or more positions. Suitable optical Teceiving devices 
would be a photodiode, an optical fibre coupled to a photodiode or a waveguide 
coupled to a fibr«. Let (x,y) be the Go-ordiaate, system at fee Fourier plane. The 
origin for x and y is the position where the lens optical axis intersects the 
Fourier plane. Let the Fourier transform of Hd(u.v) f(u,v) exp i. (p(u.v) be Fo(x,y) 
exp i d(x,y) where Fd(x,y) describes the amplitude and 6(x,y) describes the 
phase. Hence the measured intensity is proportional to the term Fo(x,y). 

Now change the hologram pattern in a known way at the position (uq, voj where 
it is required to characterise the beam phase and amplitude such that the 
hologram pattern H(u,v) becomes that shown in equation (2): 

H{u,v) = («, v) + (m.v) at (u,v) close to (uo, Vo) 

(2) 

H (m, v) = Hq («, v) elsewhere 

Hence in a known neighbourhood of the point (uo ,vq) there is a perturbation in 
the hologram pattern. Therefore there is also a perturbation in the output field 
from tiie SLM, given by Hi(u,v)f(u.v) exp i ^u,v). The iucident field amplitude 
and phase may be approximated to be uniform across the perturbation region, 
in which case the perturbation in the output field from the SLM may be 
represented as Hj(u,v) f[uo,vo) exp i (Muo,vo). As an example consider the 
perturbation to be a uniform change in phase modulation over the area of a 
single square pixel of side p, from an initial phase qo to a new phase qi. This 
perturbation may be described as a "flashiug pixel". For this example Hj(u,v) 
may be represented by equation (3) 

H,(M,v)=expigi-expi^o V(tt.v) :|«-Uo|^-|-.|v-Vo|<-| 
= 0 elsewhere 

(3) 
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Let the Fourier transfonn of Hi(u,v) be g(x,y)exp i yK^,y). Hence the Fourier 
tfansfSfm Of ni{VL,v) ftUo,Vo) exp i 4)Cuo,vo) is given by f[TidiVd) exp i 0(tto.V6) 
g(x,y) exp i y/(x,y). So this perturbation field contains information about the 
phase and amplitude of the incident beam. For the example given above, the 
perturbation field is given by equation (4). 



(4) 

exp i yrix, y) exp z^K ' ) = ^ 



where /is the focal length of the Fourier lens. From the equation the amplitude 
of the perturbation field in the Fourier plane is proportional to the amplitude of 
10 the incident field at the perturbation in the hologram. The phase of the 
perturbation field includes a constant component equal to the phase, of the 
incident field at the perturbation in the hologram, and also a linear component 
in X proportional to the position of the perturbation in tibie hologram. 

15 Returning to the general case, in the Fourier plane the field is the Fourier 
transform of the output field from the SLM. Given that a Fourier transform is a 
linear operation the field is the sum of the individual Fourier transforms of the 
original and perturbation field from the SLM. Hence the total field in the 
Fourier plane, F(x,y) is that given by (5): 

20 

F(jc, y) = Fo(jc,y)expf6^(jc, y) + /(Mo,Vo)expi(K«o>Vo)g(J^, y)expiXx,y) 
(5) 

The analysis given hereafter assumes the optical output to be detected by a 
25 single photodiode that is small enough such that the field amplitude and phase 
in the Fourier plane may be considered uniform over the active area of the 
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photodiode. Equations are derived for the response and data fitting methods are 
described to measure the amplitude variation and phase variation of the field 
incident on the. spatial light modulator. As mentioned previously, other 
receiving elements could be used. Examples are an optical fibre (single mode 
or multimode) coupled to a photodiode. In the case of a single mode fibre a 
mode stripper should also be used. Other example receivers are a larger 
photodiode or an array of photodiodes. For each case, knowing the physics of 
the receiving process (which involves a coupling efficiency calculation for the 
optical fibre case) analytical expressions, may be used to derive the receiver 
response, given the incident field as described in equation (5). As is 
demonstrated for the case of a small single photodiode, data fitting methods 
may be derived, based on said analytical expressions, to measure the amplitude 
.variation and phase variation of the field incident on the spatial light 
modulator. 

Assuming the field is detected directiy by one or more photodiodes, the 
induced photocurrent at position (x,y) is proportional to the local intensity. The 
intensity at the Fourier plane contains 3 terms. 

The expression for the intensity, I(x,y)y is given by equation (6). 
Iix,y)ccF^\x,y) 

+ 2Fo (x, y)g(x;y)f(Uo , ) cos{^(ao , v,,) + i/r^x, y) - 9{x, j;)} (6) 
+/^(«o'Vo)5"(-«,3') 

The first term , F^ix^y) is the original intensity, before the perturbation was 
applied.' . 



9 



The second tenn 2Foix,y)gix,y,)fiu„Vo)cosmuo,v^)+vrix,y)-0ix,y)}is a 
coherent coupling term between the original field in the Fourier plane, and the 
field component created tee by the perturbation at Uie SLM. IMs second term 
contains information about the phase and amplitude of the beam incident on the 
5 SLM. For the example hologram perturbation the second term also contains 
information about the flashing pixel. 

The liiird term f\uo,v„)g''ix, y) is the intensity that would appear at the Fourier 
plane if the perturbation was acting on its own, with the original field removed. 
10 This third term contains information about the amplitude of the beam incident 
on the SLM. but not the phase. The values that are to be extracted from the 
system are f(uo. vo) and </i(uo, vo). The other "unknowns" are g(x,y) and 
Wix,y)-0(.x^y}- 

15 For a general case where the incident beam phase varies in an unknown way. 

Usually what is desired is the relative amplitude variation of the beam, but not 
the absolute ampUtude. The terms Fo(x.y)g(x,y) and g(x,y) are independent of 
the position of the perturbation at the SLM, so can be considered as a 
multiplying constant. The ampUtude/fwo, vo) may be measured in several ways. 

20 

A first method is to measure the intensity in a region of the Fourier plane where 
Fo(x.y) is very weak but g(x,y) is relatively stronger such that only the third 
term is significant. The region of the appUed perturbation at the SLM is broader 
than the region occupied by the incident beam, hence at the Fourier transform 
25 plane the region occupied by the Fourier transform of the perturbation is 
narrower than the region occupied by the Fourier transform of the original 
field. Hence by varying the position on the SLM where the perturbation is 
appUed a set of measurements of (g(x.y)f(uo,vo)f may be built up at different 
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positions (uo,vo). By takmg the square root of these measurements values for 
the relative field amplitude at these positions may be derived. 

In general there may be sidelobes or spurious diffraction orders in the beam 
Fo(x,y) making it difficult to find such a region. Often it may be that g(x,y) is 
relatively weak and Fc(x,y) is stronger. Therefore it becomes appropriate to use 
the second term to measure the amplitude. In such a region the third term may 
be neglected and the second term estimated by subtracting the receiver output 
before the perturbation was appUed firom the receiver output in the presence of 
the perturbation. 

A second method is to step through a sequence of distributions Hd(x.y) chosen 
to have the same values for Fo(x,y) and known changes in the values for 6(x,y). 
In an embodiment example Ho(x.y) is a binary phase pattern, in which F(^x,y) is 
independent of the relative position of the pattern on the SLM, but e(x.y) 
changes in a known way as the pattern position is changed on the SLM. In 
general let 6(x,y) be expressed as given by equation (7). 

(7) 

where m is a variable that represents the known phase associated with the 
pattern position, and 0o represents the sum of the background phase component 
from ithe hologram tiiat is not affected by the pattern position and the phase 
from the incident field. 

By applying a set of patterns, at different relative positions on the SLM a set of 
values may be built up for the term ^ven in equation (8): 
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2Fo(;c. y)g(x, y)/(«o. Vo) cos{(*(Uo ,Vo) + V^ix, y) - O^ix, y) - 0(m)} 
(8) 

which may be further expressed as given by equation (9). 

cxix, y) nuo,vo)cos{j3(x, y,u„,Vo)}cos{e(pi))+aix, y)/(Uo.Vo)sin{^(x, y,Mo,Vo)}sin(6'(m)) 
(9) 

where c«fxiy; = 2Fo(x,y)g(x,y) and^x,y,Uo,vo) = (p(uo.vo)+V^x,y)-eo(x,y). 

Equation (9) may be considered as a linear equation in cos(6(m)) and sin(d(m)), 
with unknown coefficients c = cx(x.y) cos /3(x,y,uo,vo) f[uo,vo) and d = cx(x,y) sin 
/3(x,y,uoyo) f(uo.vo). Any suitable data fitting method may be used to extiract 
values for these coefficients c and d. A relative value for the amplitudd-may 
then be calculated from equation (10) (remembering tiiat a(x,y) acts -like a 
multiplying constant at any Fourier output position (x^y)). 
«(A:,3')/("o.Vo) = Vc^ + d' 

(10) 

Other methods could also be used to process the data. 

While the above exaiiq)le assumes an invariant Fd[x.y) and known changes in 
e(x,y) hologram pattems could be selected that maintain an invariant e(x,y) and 
change Fd(x,y)m a known way, or that change both in a known way. Another 
method is to change the amplitude and/or the phase of the perturbation 
hologram in a known way. The general principle is to be able to take a 
sequence of measurements with some known parameters changing in a known 
way and use data fitting to extract the values of the unknown parameters. These 
unknown parameters may be combined in a way that extracts the incident field 
ampUmde without requiring knowledge of the incident field phase. 
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Another mfiasurement diat inay be required is to check that the pixels are 
functioning correctly. Take a measurement of the output intensity with the 
original hologram applied: this could be the coupling efficiency into a fibre, or 
the output from a photodiode, or the outputs from a photodiode array. Now 
flash a pixel, one at a time, preferably with a phase perturbation qi - go close to 
Tc so as to maximise the effect. Calculate the relative amphtude as described 
above. If there is no significant ampKmde then either the pixel is outside die 
area of the incident beam or the pixel drive circuit is not functioning. If the 
ampUtude varies in a smooth manner hke the expected profile then the pixels 
are working correctly. If there is a discontinuous jump in the response to a 
higher value then the pixel may be shorted to another pixel. If there is a 

discontinuous jump to a lower value then the pixel is not being driven 

. correctly. 

In some systems what is required is to check or measure the position of the 
field with respect to the SLM pixels. The peak measured amplitude wiU occur 
when the pixel being flashed coincides with the peak of the incident beam. 

In many cases it may also be required to measure the phase variation in the 
incident field. The methods described hereafter explain how to measure the 
variation in the u direction in the interests of clarity. However the method is 
easily extended to measurements in both u and v directions. 

Given the data fitted values of the coefficients c and d the phase tenn 
fi(x,y,uo.Vo) may be calculated using equation (1 1). 

COs(;g(^,y.U„.Vo))= sl( 7 

sm{/3{x,y.u„v,))=-. ^ 

(11) 
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Let the initial position (uo, vq) be a reference point such that the phase at any 
arbitrary point (u,v) is measured with respect to the phase a.t(uo, vq). Using (4) 
there is obtained a genial ej^ession for ^x,y,uo,vo) as shown in equation 
(12): 

5 iS{x, y.Ko .vo) = — ^ + ?>(«o . vo )+ 2" ° J 

(12) 

Now consider the effect of changing the position of the hologram perturbation 

to (u,Vo) but maintaining the position of the output receiver at (x,y). The value 

of P.becomes that shown in equation (13). 

of \ 2^0:^ ,/ . i2m>oX . ^ . gp+gi p /_ ^ J' 
10 fi{x,y,u,Vo)=-^+^{u,Vo)+^^-^^+—+— <9o(^.3')j 

(13) 

The difference in the values of p contains two terms, the first (in; curly 
brackets) is the required phase difference being measured, while the secolid is a 
linear phase term, as shown in equation (14). 

15 ^{x, y, M. vo )- J3{x, y.uo . vo ) = {^(m, Vp )- <»(ko . Vq )}+ ^ "' 

(14) . 

Therefore to extract the required phase difference a method is needed to 

1' 

estimate or measure the linear phase tenn. Given knowledge of the pixel pitch 
on the SLM it is straightforward to calculate (u-uo). Similarly knowledge of the 
20 wavelength, A,, is established. The focal length can be measured. Therefore it is 
possible to know to a high degree of accuracy the value of 2:/i(u-uo)/f/i What 
may be difficult is to obtain an accurate value of x, the distance in the x 
direction fi-om the point where the lens optical asds intersects the Fourier plane 
to the receiving element. 

25 
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One method to overcome this difficulty is to pre-calilnrate the system with a 
reference beam. This could be a well-collimated beam, or a beam of known 
phase variation (-using for example a- pioliole to generate the reference beam), 
or a beam such that the phase of any arbitrary beam needs to be measured with 
respect to the reference beam. Let the subscript R represent the corresponding 
measurements for tihe reference beam. Hence equation (15) is obtained: 

yg.(^.y.".vJ-^^(;c.3>.«o.Vo) = kkvJ-(^>,(«„.v„)}-H ^^"7<')^ 

(15) 

Therefore if the p difference is measured for the reference beam and the phase 
variation across the reference beam is known, the value of the linear phase term 
may be calculated. Furthermore, by repeating this process for several different 
positions u on the SUM a data fitting method may be used to calculate the value 
of the unknown parameter x. 

Without knowledge of the phase variation of the reference beam but with a 

need to use it as a baseline, it is necessary to measure the >9 difference for the 

beam under test, and subtract from it the >^ difference for the reference beam, at 

the same positions u and Uq. Alternatively interpolation may be used to predict 

the ^difference for lie reference beam if it has not been measured at precisely 

the required points. An expression may be obtained for the phase variation with 

respect to the baseline variation, as shown in equation (16): 

^(jc,3'.«.Vo)-^(j:,>',Mo.Vo)-{/^/e(^,}',«!vo)->ff«(x,3;,«o,Vo)} 
= Vo)-(^(«o , Vo (m, Vo )- j*^ ,Vo )} 

(16) 

It is possible to perform the invention with embodiments using either 
multiphase or binary phase SLMs. Multiphase pixellated SLMs may 
incorporate an integral or non-integral quarter-wave plate or a wave plate 
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having a similar effect to provide polarisation insensitivity where a liquid 
crystal having out of plane tilt is used. 

The above description refers to LCOS SLMs. However the invention is not so 
limited but instead extends to the full scope of the appended claims. 
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1. A method of measuring amplitude and phase variations in a beam of 
light comprising causing the beam to be incident upon a spatial airay 
displa3dng a pixellated first phase distribution, in a measuring region of said 
spatial array causing the phase distribution to assume a new value while 
retaining the first phase distribution outside the measuring region, in the 
Fourier plane determining the change in intensity resulting from the change in 
phase distribution 

2. A method of characterising a beam of light, comprising 
disposing a LCOS SLM in the path of the beam; 
causing the LCOS SLM to display a first hologram pattern; 

at a location in said beam where the ainplitude and phase of the beam 
are to be characterised, changing the hologram pattern to a second hologram 
pattern; and 

measuring the effect of said change by measuring an intensity. 

3. A method as claimed in claim 2, wherein, the output jfrom the SLM is 
measured in the Fourier plane to detect the Fourier output. 

4. A method as claimed in any preceding claim, comprising measuring the 
intensity in a region of the Fourier plane where the amplitude distribution 
associated with the beam modulated by the original hologram is relatively 
stronger, but the amplitude distribution at the Fourier plane of the field 
component created by the perturbation in the hologram is relatively stronger, 
varying the position on the SLM where the perturbation is applied. 
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5. A method as claimed in claim 4 further comprising taking the square 
root of a set of values obtained. 

6. A method as claimed in any preceding claim, comprising stepping 
5 through a sequence of phase distributions. 

7. A method as claimed in any preceding claim, comprising varying the 
phase shift in a respective single pixel. 

10 8. A method as claimed in claim 3 comprising manipulating hologram 
patterns to obtain information related to a coherent coupling term to thereby 
derive amplitude information. 

9. A method as claimed in claim 3 comprising manipulating hologram 
15 patterns to obtain informatibn related to a coherent coupling term to thtereby 

derive phase information. 

10. Apparatus for measuring ampUtude and phase variations in a beam of 
light, the apparatus comprising a pixeUated spatial array, each pixel being 

20 controllable to apply any of plural phase shifts to input light, whereby the array 
displays a desired distribution of phase modulation, 

means for causing the array to display a first selected distribution of 
phase modulation; means for changing the first distribution in a measuring 
region of said spatial array to assume a new distribution while retaining the 
25 first phase distribution outside the meas\iring region, 

means disposed in the Fourier plane for determining a change in 
intensity of light resulting from the change in phase distribution 
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1 1. Apparatus as claimed in claim 10, wherein the spatial array has only two 
possible values of phase shift per pixel. 

12. Apparatus as claimed in claim 10, wherein the spatial array has more 
than two possible values of phase shift per pixel. 

13. Apparatus for characterising a beam of light, comprising 

a LCOS SIM arranged so that a said beam of Ught can be incident upon 

it; 

means for causing the LCOS SLM to display a first hologram pattern; 
means for changing the hologram pattern to a second hologram pattern at a 
location in said beam where the amphtude and phase of the beam are to be 
characterised; and 

means for measuring an intensity of Ught to detennine Hie effect of said 
change of hologram pattern. 

14. Apparatus as claimed in claim 13, wherein the means for measuring is 
disposed in the Fourier plane to detect the Fourier output 

15. Apparatus as claimed in claim 13 or 14, further comprising a lens for 
providing the Fourier output 

16. Apparatus as claimed in claim 13, 14 or 15, further comprising a mirror 
for providing the Fourier output 
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